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1 . INTRODUCTION 

A model  of  high  frequency  radio  (HF)  skywave  modes  under  ambient 
conditions  is  basic  to  the  study  of  HF  system  response  to  disturbed  condi- 
tions, such  as  those  which  would  be  encountered  in  a nuclear  conflict. 

Traditionally,  HF  skywave  vulnerability  has  been  calculated  under 
the  assumptions  that  only  D-layer  absorption  is  important,  that  ambient 
propagating  modes  are  an  adequate  approximation  to  modes  which  propagate 
after  nuclear  explosions,  and  that  a careful  description  of  late  time  (ten 
minutes  to  hours)  debris  motion  is  unnecessary  so  long  as  an  assured  defense 
conservative  model  is  used. 

Recent  interest  in  adaptive  HF  systems  for  the  Air  Force  and 
MEECN,  as  well  as  interest  in  the  closely  related  class  of  meteor  burst 
communication  systems,  has  forced  a re-evaluation  of  all  three  of  the  above 
assumptions.  More  realistic  models  of  late  time  nuclear  debris  motion  (the 
source  of  nuclear  radiations  which  ultimately  cause  D-layer  absorption) 
have  recently  been  developed1’2  to  allow  improved  absorption  calculations. 

In  this  report  we  describe  a method  of  ambient  HF  skywave  mode  calculation 
designed  to  readily  accept  disturbed  ionospheric  conditions,  and  thus  to 
allow  inclusion  of  some  nuclear  disturbances  important  to  propagation 
geometry,  necessary  for  more  accurate  calculation  of  D-layer  absorption 
effects  as  well  as  for  removal  of  the  second  and  third  traditional  assump- 
tions quoted  above. 
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The  necessity  for  development  of  an  ambient  skywave  propagation 
mode  calculation  requires  justification.  Elaborate  techniques  for  detailed 
calculation  exist  in  the  scientific  community,3  which  have  been  adapted 
into  a nuclear  effects  code.1*  Less  elaborate  techniques  exist  for  long 
term  average  mode  prediction,5  to  which  simplified  fits  have  been  made, 
and  these  fits  are  commonly  used  in  the  nuclear  effects  community.6  To 
understand  the  requirement  for  a technique  of  intermediate  complexity  re- 
quires a brief  discussion  of  anticipated  nuclear  effects  of  HF  propaga- 
tion in  relation  to  characteristics  of  available  calculational  techniques. 

Anticipated  nuclear  effects  which  may  conveniently  be  taken 
into  account  by  alteration  of  normal  ionospheric  parameters  are  changes 
in  ambient  layer  (i.e.,  F2,  FI,  or  E-layer)  electron  density,  thickness, 
height,  and  tilt  due  to:  (1)  propagating  acoustic  gravity  waves,  (2)  ducted 
waves,  (3)  burst  region  disturbance  (F-layer  depletion,  etc.),  or  (4)  bomb 
or  debris  irradiation.  Another  set  of  effects  exists  (M-modes,  fireball 
modes,  plume  modes,  etc.)  which  are  more  conveniently  calculated  indepen- 
dently and  which  will  not  concern  us  here.  To  calculate,  even  crudely, 
effects  of  the  above  list  on  propagation  geometry,  thus  on  existence  and 
efficacy  of  ambient-like  modes  and  creation  of  new  modes,  the  ambient 
mode  calculation  must  readily  take  into  account  ionospheric  variations  all 
along  the  propagation  path  in  terms  of  variability  of  reflection  height 
and  ionospheric  tilt,  both  along  and  across  the  great  circle  path.  In 
fact,  non  great  circle  propagation  must  be  treated.  The  current  ITS78 
code,5  the  community  standard,  evaluates  ionospheric  conditions  only  at  the 
great  circle  path  center,  derives  from  that  a single  reflection  height, 
and  does  not  consider  tilt  or  off- great  circle  propagation.  This  is  satis- 
factory for  long  term  average  predictions,  the  objective  of  that  code,  but 
does  not  lend  itself  to  situations  similar  to  nuclear  effects  where  a portion 
of  the  path  is  strongly  disturbed.  The  WRECS  code,6  using  a fit  to  ITS78 
predictions,  allows  for  two  values  of  reflection  height  corresponding  to 


4 


day  or  night  but  is  otherwise  restricted  to  ITS78  predictions.  The  NUCOM 
code4  uses  ray  tracing  techniques  and  thus  can  be  much  more  versatile,  but 
such  techniques  are  expensive  for  our  purposes. 


The  method  developed  here  to  meet  requirements  of  nuclear  effects 
studies  is  approximate  but  more  accurate  than  justified  by  predictability 
of  either  the  ambient  or  nuclear  disturbed  ionosphere.  It  is  more  expen- 
sive in  computer  time  than  one  would  prefer  but  still  vastly  cheaper  than 
normal  ray  trace  techniques;  sophisticated  in  that  it  accounts  for  ionos- 
pheric variation,  tilt,  and  off-great  circle  propagation;  but  in  its  pres- 
ent implementation  rudimentary  in  that  only  two  modes  can  be  calculated  for 
a given  number  of  hops,  no  extraordinary  ray  modes  are  calculated,  etc.  The 
restriction  to  two  modes  per  hop  number  could  be  readily  remedied  and  this, 
at  least,  is  hoped  to  be  done  soon.  The  basic  technique  is  believed  to  be 
sufficiently  versatile  that  extension  toward  more  realistic  treatment  is 
feasible.  As  reported  here  the  method  has  been  develoepd  to  a useful  state, 
one  that  is  more  realistic  than  others  in  many  important  respects,  as  well  as 
one  that  is  reasonably  inexpensive  in  application. 

An  outline  of  the  overall  calculation  is  given  in  the  next  sec- 
tion and  specific  components  discussed  in  following  sections. 
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2. 


OUTLINE  OF  THE  CALCULATION 


The  calculation  is  naturally  broken  into  two  logical  modules, 
each  of  which  consists  of  several  component  subroutines. 

2.1  Ionospheric  Data 

The  first  module  is  exercised  once  only,  before  any  propagation 
calculations  are  attempted.  It  uses  an  ionospheric  model  to  create  a data 
table  of  ionospheric  parameters  for  each  desired  transmitter-receiver  link. 
Each  such  data  table  consists  of  three  sets  (one  for  each  ionospheric  layer, 
E,  FI,  F2)  of  three  parameters  describing  parabolic  ionospheric  layers 
above  several  points  on  the  great  circle  path  at  a number  of  times.  The 
three  layer  parameters  are:  (1)  maximum  electron  density,  (2)  thickness 
of  the  layer,  and  (3)  altitude  of  the  maximum  electron  density.  The 
spatial  points  are  arbitrary  in  number  and  spaced  equally  along  the  great 
circle  path.  The  number  of  times  is  also  arbitrary  and  determines  a number 
of  equal  segments  of  a twenty-four  hour  period  at  which  each  layer  will  be 
specified  at  each  point  along  the  great  circle.  Normally  we  choose  twenty- 
four  times,  one  for  each  hour  of  the  day  (plus  one  repeated  entry  to  allow 
the  table  to  begin  at  0 hours  and  end  at  24  hours),  and  seven  space  points; 
which  yield  six  equally  spaced  points  along  the  path  (both  end  points  are 
included),  normally  these  great  circle  points  are  separated  by  less  than 
15°  of  arc  each  (that  is,  corresponding  to  less  than  an  hour  apart  local 
time).  With  ionospheric  models  at  our  disposal 5,7  one  hour  equivalent  seg- 
ments have  been  found  to  be  about  twice  as  dense  as  necessary  to  allow 
interpolation  in  both  space  and  time  to  the  limit  of  model  fidelity. 

The  purpose  of  this  array  of  (3  layers)  * (3  parameters)  x 
(typically  7 space  points)  x (typically  25  time  points)  = 1575  numbers 
per  link  is  to  allow  rapid  interpolation  in  space  and  time  for  ionospheric 
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quantities  without  having  to  repeatedly  calculate  potent!  ->lly  complex  ionos- 
pheric models  for  identical  or  similar  conditions.  Heights  and  thicknesses 
are  interpolated  linearly  between  adjacent  points  and  electron  density 
logarithmically.  In  Section  5 we  shall  show  that  despite  appearances, 
this  form  of  ionospheric  data  table  does  not  restrict  consideration  to 
great  circle  paths. 

A typical  comparison  of  the  three  parabolic  layer  ionospheric 
approximation  used  here  to  the  ITS78  model,5  the  Aerospace  Model7  and  to  data7 
is  shown  in  Figure  1.  Cases  can  be  found  in  which  the  three-layer  approxima- 
tion fits  the  model  somewhat  better  and  somewhat  worse  than  Figure  1.  It  is 
also  easy  to  find  both  better  and  worse  comparisons  between  the  models  and 
data.  Figure  1 illustrates  the  general  result,  which  appears  to  be  valid 
for  all  ionosphere  models  which  might  be  used  as  a basis  for  the  fit,  that 
error  resulting  from  the  three  parabolic  layer  approximation  is  small  com- 
pared to  error  due  to  the  uncertain  behavior  of  the  real  ionosphere. 

2.2  Mode  Calculation 

The  second  module  calculates  mode  geometry.  Inputs  to  this  module 
are  a specified  time,  ray  type,  link,  frequency,  and  number  of  hops.  All  of 
these  quantities  are  arbitrary  except  for  obvious  restrictions,  frequency 
and  number  of  hops  must  be  greater  than  zero;  and  except  for  a bow  to  the 
nature  of  Fortran  compilers,  presently  we  dimension  for  no  more  than  eight 
hops.  The  ray  type  specifies  either  the  low  ray  or  high  (Pederson)  ray. 
Outputs  of  this  module  are  either  specification  of  mode  geometry  or  infor- 
mation concerning  the  reason  no  such  mode  exists. 

Mode  geometry  is  specified  in  terms  of  the  ends  of  a number  of 
straight  line  segments  (passes)  connecting  reflection  points.  A one-hop 
mode  is  specified  by  latitude,  longitude,  and  altitude  above  the  earth's 
surface  of  three  line  segment  ends;  (1)  the  transmitter  coordinates. 
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Altitude,  km 


Figure  1.  Typical  ambient  ionosphere.  The  HFNET  three  parabolic  layer  fit 
is  compared  to  the  Aerospace  model  from  which  it  is  derived,  to 
the  ITS78  model,  and  to  data  for  1200  local  time  above  Arecibo, 
Puerto  Rico,  Dec.  1968.  Smoothed  Zurich  sunspot  number  is  110. 


(2)  the  equivalent  triangular  ionospheric  reflection  point,  and  (3)  the 
receiver  coordinates.  A two-hop  mode  requires  specification  of  coordi- 
nates for  five  points;  the  two  end  points,  two  ionospheric  reflection 
points,  and  an  intermediate  surface  reflection  point.  In  general  three- 
space  coordinates  must  be  specified  for  2n  + 1 points  for  an  n-hop  mode. 

A mode  fails  to  exist  if:  (1)  the  F2  layer  fails  to  reflect 
the  ray  on  any  hop,  (2)  transmitted  launch  angle  for  the  mode  is  below 
a specified  value,  (3)  received  angle  is  below  a specified  value, 

(4)  ionospheric  tilt  causes  the  ray  to  be  reflected  back  upon  itself,  or 

(5)  ionospheric  tilt  causes  the  ray  to  be  deflected  forward  at  too  great 
an  angle. 

Steps  in  calculation  of  n-hop  mode  geometry  are  as  follows: 

1.  The  great  circle  path  is  divided  into  n equal  segments, 
one  for  each  hop. 

2.  For  each  segment  a one-hop  mode  is  calculated  according 
to  local  ionospheric  conditions  as  follows: 

a.  Assume  reflection  occurs  at  the  F2-layer;  calculate 
the  equivalent  triangular  reflection  altitude,  hr> 
for  a parabolic  layer  by  iteration  patterned  after 
Reference  5. 

b.  If  successful,  calculate  refraction  effect  due  to 
the  FI  layer.  If  the  FI  layer  reflects  the  ray, 
carry  out  (a)  for  the  FI  layer. 

c.  Calculate  refraction  effect  due  to  the  E-layer.  If 
the  E- layer  reflects  the  ray,  carry  out  (a)  for 
the  E- layer. 
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The  situation  at  the  end  of  step  (2)  is  depicted  in  the  sketch, 
provided  all  hops  were  successfully  reflected.  At  this  point  in  the  cal- 
culation the  mode  is  specified  by  latitude  and  longitude  end  points  equally 
spaced  along  the  great  circle  path  but  each  ionospheric  reflection  height 
is  different,  a violation  of  the  laws  of  optics. 

3.  Adjust  the  pass  ranges,  r^,  so  that  incident  and  reflected 
angles  are  equal,  accounting  for  differences  in  h^. 

4.  Calculate  ionospheric  tilt  angles  both  along  the  great 
circle  and  across  the  great  circle  at  each  ionospheric 
reflection  point. 

5.  Adjust  all  pass  coordinates  (latitude,  longitude,  altitude) 
to  account  for  tilt  along  the  great  circle  to  first  order 
in  tilt  angle  (that  is,  in  the  vicinity  of  each  reflection 
point,  the  tilt  is  constant).  If  the  ray  is  reflected  back 
upon  itself  at  any  point,  or  if  the  cumulative  effect  of 
all  tilts  cause  the  ray  to  land  too  far  beyond  the  re- 
ceiver, reject  the  mode. 
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6.  Calculate  off-great  circle  displacements  of  all  pass  coordi- 
nates and  new  reflection  altitudes  due  to  across  great  circle 
tilt . 

7.  Find  launch  and  received  ray  elevation  relative  to  the  horizon, 
and  azimuth  relative  to  the  great  circle  path.  If  either 
elevation  is  below  a specified  limit,  reject  the  mode. 

At  the  end  of  this  process  an  approximation  has  been  found  to  a 
"good"  propagating  mode  with  specified  number  of  hops,  within  the  limitations 
indicated  above. 

Possible  paths  which  include  hops  of  the  form  indicated  in  the 
sketch  are  not  found.  The  situation  depicted  in  (a)  is  rejected  because 


r2  is  negative,  the  ray  was  reflected  back  upon  itself.  Within  the  mathe- 


matical restrictions  of  first  order  tilt  angle  treatment  this  situation  can 
only  arise  when  hr  is  greater  than  the  chord  between  the  two  ground  points, 
and  can  only  represent  a reasonably  plausible  situation  when  the  second  point 
represents  receiver  or  transmitter.  These  conditions  are  too  severe  for  am- 
bient skywave  propagation.  Situation  (b)  is  rejected  because  the  calculation 
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does  not  allow  for  existence  of  the  second  reflection  point,  2,  and  assumes 
that  the  ray  continues  along  the  dotted  path.  This  situation  seems  unlikely 
under  ambient  conditions  but  could  occur  in  a disturbed  environment.  It 
must  be  handled,  if  at  all,  in  future  as  a "disturbed"  mode,  with  appro- 
priate logic  specified. 

The  reader  will  note  that  no  step  between  (1)  and  (7)  was  iterated. 
The  only  iterations  occur  in  finding  single-hop  solutions  to  reflection  from 
a single  parabolic  layer  in  step  (2).  With  this  exception  the  mathematical 
treatment  runs  "open  loop"  and  could  be  made  more  rigorous  by  iterating  on 
each  refraction  correction,  after  range  correction  for  differing  reflection 
heights,  and  after  corrections  for  tilt.  Such  a procedure  would  be  pro- 
hibitively expensive  in  calculation  time  and  would  produce  numerical  results 
significantly  different  from  the  present  technique  only  for  cases  where  some 
reflection  occurred  very  near  a critical  point  (e.g.,  almost  or  barely  re- 
flected off  some  layer).  Predictability  of  either  natural  or  disturbed 
ionospheric  conditions  is  not  sufficiently  accurate  to  justify  such  fine 
tuning. 


Of  much  greater  importance  is  the  fact  that  this  method  does  not 
calculate  all  possible  "sensible"  rays.  This  defect  is  illustrated  in  the 
sketch,  in  plane  geometry  for  illustrative  purposes. 
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Our  technique  first  finds  a tentative  (low  or  Pederson)  ray  which 
reflects  off  the  highest  available  layer,  the  dashed  ray  of  the  sketch  which 
reflects  in  the  vicinity  of  point  A.  Next  the  refraction  angle  due  to  pas- 
sage through  the  lower  layer  at  point  B is  calculated.  If  refraction  is 
not  large  enough  to  reflect  the  ray  at  B,  a correction  to  the  effective 
triangular  reflection  point  A is  introduced  and  calculation  proceeds  to  the 
next  step.  At  this  point  our  technique  has  irretrievably  missed  a potential 
lower  ray  which  reflects  at  C.  Under  almost  all  ambient  conditions  the  lowest 
ray  suffers  greater  absorption  and  also  is  a more  likely  candidate  for  rejec- 
tion due  to  low  launch  angle,  thus  is  not  likely  to  be  useful  or  to  cause 
multipath  problems.  It  is  unusual  but  possible  that  nuclear  effects  could 
preferentially  interfere  with  the  upper  ray  enough  to  make  the  ray  at  C the 
best  mode.  This  is  a weakness  in  current  logic  which  probably  should  be  fixed. 

On  the  other  hand,  if  the  ray  is  reflected  at  B then  the  tech- 
nique shifts  attention  to  the  lower  layer  and  calculates  only  the  ray 
which  reflects  at  C.  Another  option  would  be  to  search  at  higher  angles 
for  a ray  which  is  strongly  refracted  at  D,  closer  to  the  critical  angle 
for  the  frequency  if  such  exists,  and  which  traverses  the  path  indicated 
by  the  light  solid  line,  a high  ray. 

Thus  the  technique  will  select  preferentially  the  dashed  ray,  then 
the  heavy  ray;  but  not  both. 

The  only  time  the  method  selects  a ray  like  that  refracting  at 
D is  in  the  rare  instance  when  point  B happens  to  lie  very  near  the  critical 
angle  for  the  lower  layer.  This  feature  is  generally  good  because  when 
refraction  occurs  in  the  vicinity  of  the  critical  angle  a significant  de- 
focusing  occurs  with  attendant  loss  of  signal  strength.  But  if  strong 
D-layer  absorption  occurs  on  the  heavy  ray,  the  defocusing  loss  at  D could 
be  less  than  absorption  loss;  in  this  event  the  method  misses  the  best 
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propagating  path.  Again,  a more  elaborate  procedure  could  be  devised  to 
cure  this  weakness. 

Other  cases  of  interest  which  are  not  handled  at  present  include 
spread  F and  sporatic  E multipath,  and  cases  where  unusual  spatial  depen- 
dence of  the  layers  produce  freak  (but  occasionally  real)  propagation  paths. 

In  summary,  the  selected  technique  is  believed  to  produce  the 
strongest  propagating  mode  under  most  conditions,  but  has  weaknesses  in 
that  occasionally  a potentially  significant  alternate  mode  is  missed. 

The  most  important  improvement  would  be  to  include  sporatic  E and  spread 
F effects.  These  can  be  treated  to  an  adequate  approximation  for  many 
cases  as  a perturbation  on  the  successful  ray  without  incurring  a signi- 
ficant penalty  in  computation  time  or  logic.  Treatment  of  other  alternate 
rays  will  incur  a noticeable  penalty  in  computation  time,  and  to  some 
extent  in  complexity  of  logic. 


3. 


SINGLE  HOP  CALCULATION 


After  the  great  circle  path  for  an  n-hop  mode  has  been  divided 
into  n equal  segments,  each  hop  is  calculated  independently,  as  indicated 
in  steps  (1)  and  (2)  of  Section  2.  Here  we  discuss  step  (2),  calculation 
of  a single  reflection  from  a single  parabolic  layer,  and  correction  for 
refraction  in  passage  through  a single  parabolic  layer. 

3.1  Reflection 

Our  technique  differs  from 
that  of  Reference  5 substantially 
only  in  rearrangement  of  their  iter- 
ation scheme.  For  more  details  the 
interested  reader  is  referred  to 
Reference  5. 

The  scheme  is  a prediction- 
correction  method  based  on  the  quan- 
tity 

x = occ  cos0  , 

where  a is  the  ratio  of  broadcast  frequency  to  critical  frequency  at 
reflecting  layer  maximum,  The  factor  c is  a correction  for 

ionospheric  curvature,  and  0 is  defined  in  the  sketch. 

Let  hQ  be  altitude  of  the  layer  bottom  and  a be  layer 
thickness,  then  the  apparent  reflection  height  ("equivalent"  or  "equivalent 
triangular"  height),  corresponding  to  an  assumed  value  of  x is 


15 


h 


(2) 


To  start  the  iteration  we  calculate  three  values  of  x given 
initial  guesses 

x(xx)  , x(x2)  , x(x3) 

Assuming  the  calculated  value  of  x to  be  a quadratic  function  of  the  input 
value,  we  can  solve  for  two  values  of  input  x which  should  reproduce  them- 
selves, 
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x = x ± /F  /B 
n o v o 


where 

Xo  = I tXl  ‘ x2  " R(x3  ' x2^/[(xi  " x2^  " R(x3  ' x2^ 
R = [x(x1)  - x(x2)]/[x(x3)  - x(x2)] 

B = [xfXj)  - x(x2)]/[x^  - x2  - 2xq(x1  - x2)] 


Fq  = - B(X1  - Xo^ 


If  F > x then  no  solution  exists,  the  ray  penetrates  the  ionosphere, 
o o 


It  is  necessary  to  place  a logical  guard  against  prediction  of 
negative  x;  if  this  occurs  set 

x = x /2  . (8 

n n-1 

It  is  also  necessary  to  guard  against  prediction  of  x^>  1,  if 
this  occurs  set 


*„  ■ (1  - Vi>/2 


We  use  x:  = 0.1,  x2  = 0.5,  x3  = 0.9  and  require  for  convergence 


|x(xn)  - xnl  < 10 


So  far,  no  case  tried  has  required  as  many  as  ten  iterations  to 


converge . 
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Refraction  by  Passage  Through  a Lower  Layer 


This  correction  is  also  based 
upon  Reference  5.  We  seek  the  change  in 
launch  angle.  A,  caused  by  the  presence 
of  a parabolic  layer  of  thickness  a^ 
whose  maximum  electron  density  occurs  at 
altitude  h. . The  previous  solution  for 

Li 

h"  allows  determination  of  the  unper- 
turbed launch  angle,  g.  Then,  from 
the  sketch  the  crossing  angle  is  given  by 


siny  = cosg/ (1+h  / R ) 

Li  6 


Define 


~ ojcosy 

where  w.  is  the  maximum  critical  frequency  of  the  refracting  layer  and 

Li 

u is  broadcast  frequency. 

If  u £ 1 the  lower  layer  reflects  the  ray  and  one  must  re- 
calculate h' . 

If  u<  1 then,  from  Reference  5 


A = 2 


aL  tanY  f 1 0 /l+u\  .1 

+ hL)  |_2u  1111  (l-u) 


The  corrected  launch  angle  is 


r = e ♦ a 
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and  can  now  be  recalculated,  including  refraction  effects. 

This  procedure  is  repeated  for  each  layer  below  the  reflecting 

layer. 


4.  RANGE  ADJUSTMENT  FOR  DIFFERING  REFLECTION  HEIGHTS 

The  method  chosen  is  sufficiently  accurate  for  our  purposes  but 
is  more  cumbersome  than  one  would  prefer. 


Method 


For  orientation  we  first  work  out  the  case  of  a flat  earth. 


Given  the  number  of  hops  and  an  arbitrary  set  of  reflection 
altitudes,  h^,  we  seek  the  launch  angle,  3,  and  the  set  of  pass  ranges, 
r^,  which  allow  3 to  be  constant  at  each  reflection  point. 

We  have 


tan3  = h^/r^ 


r.  = h./tan3 
i i 


and  the  total  range,  R,  is 


20 


So  the  launch  angle  can  be  found  from 


tanB  = R 1 h. 


(17) 


and  individual  pass  ranges  from  (15). 

With  the  solution  to  this 
easy  problem  in  mind,  we  now  proceed 
to  search  for  an  approximate  solution 
to  the  curved  earth  problem.  In  this 
case  all  of  the  surface  reflection 
angles  will  again  have  a common  value 
8,  as  illustrated  in  the  sketch  on  p 20_ 

Failing  to  find  a neat  closed  form 
solution  to  the  problem  we  proceed  to 
find  an  approximation  which  resembles 
(14)  and  is  accurate  for  h^/R«l.  Such  an  approximation  is 


h. 

tan8  = — 


[1  - rf/ (2h. R )] 
i x e 

r.  (1  + h./R  ) 
i i eJ 


(18) 


Equation  18  is  good  to  about  10  percent  for  Ik  of  3000  km  and 
r^  small  compared  to  Rg,  and  becomes  better  rapidly  as  Ik  approaches 
realistic  reflection  altitudes  of  a few  hundred  kilometers. 

Proceeding  as  in  the  plane  case,  we  solve  (18)  for  r. 


r . = 
i 


Re  ^2x^  + (1+x.)2  tan28  - (l+x^)tan& 


(19) 


where  x.  = h./R  . 

i i e 
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At  this  point  it  is  clear  that  a problem  exists,  for  if  one  sums 
Equation  19  to  force  agreement  with  the  total  range,  analogous  to  the  method 
by  which  (16)  was  obtained,  we  shall  not  be  able  to  readily  solve  for  tang. 

Accordingly,  we  make  further  approximations  by  studying  large 
6 and  small  3 limits. 


If  8 is  large  enough  (19)  can  be  expanded  to  obtain 

hi  2 

ri  (l+hVRe)tanB  + 

Then,  summing  as  before 

1 h' 
tang  R L (i+h./R  ) 
i i e 

which  looks  much  like  the  plane  case  solution.  Equation  17. 


(20) 


(21) 


This  large  B limit  occurs  in  practice  but  unfortunately  will 
not  be  exclusive,  or  even  common,  in  most  situations.  The  small  8 limit 


2R  h.  - (R  +h. ) tang  + 
e i e iJ 


(R  +h.)2  tan2g 
e i 

2\/2R  h. 

” e i 


Summing  (22)  to  find  total  range  and  solving  for  tang  yields 


(23) 
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where 


studied 

limit. 


where 


and 


»=  E V^v7 


B = E (R  »h.)2AJ2R  h. 

. e 11  e 1 


c=  E (Re+hi) 

i 


Behavior  of  the  two  limiting  solutions  can  be  conveniently 
by  choosing  a case  where  all  h^  are  equal,  then  the  large  g 
Equation  21,  becomes 


(24) 


tang  = 


Nh 


R(l+h/Re) 


(25) 


is  twice  the  number  of  hops. 


The  small  g limit.  Equation  23  becomes 


tang  = 


V2Fh 

(Re+h) 


1 - 


2R 


N\2Rh 


Define 


= yj 2N2  hRe/R2 


(26) 


(27) 


F = N(Re+h)  tang/R 
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then  Equation  25  becomes 


Fl  = q2/2  (28) 

and  Equation  26  becomes 

Fs  = q (1  ->/2/q-l)  (29) 

The  behavior  of  F vs.  q is  illustrated  in  Figure  2,  the  exact 
solution  illustrated  for  this  special  case  is  obtained  readily  from  Equation  21. 

A simple  technique  with  accuracy  adequate  for  our  purposes, 

suggested  by  Figure  2,  is  to  calculate  both  limits,  that  is,  Equations  21 

and  23,  then  choose  the  smaller  of  the  two. 

For  realistic  cases,  q is  never  as  small  as  0.1  (which  implies 
h < 1 km)  so  the  error  near  q = 0 is  not  significant.  The  maximum 
fractional  error  of  about  30  percent  occurs  at  q = 2,  typically  at  the 
very  high  launch  angle  of  30°. 

Much  of  the  error  is  removed  by  a final  adjustment  which  forces 
the  r^  to  precisely  sum  to  the  correct  range. 

4.2  Example 


An  example  which  stresses  the  algorithm  somewhat  is  summarized 
in  Table  1.  The  transmitter  is  located  at  Washington,  D.C.  and  the  re- 
ceiver at  Tel  Aviv,  at  84.92°  great  circle  distance  (9432  km).  Ionospheric 
conditions  are  appropriate  for  0600  UT,  which  places  the  western  two-thirds 
of  the  link  in  darkness  (local  time  in  Washington  is  about  1:00  a.m.  and 
in  Tel  Aviv  is  about  8:30  a.m.) 
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1 2 3 


Figure  2.  Reflection  height  adjustment.  Small  and  large  B limiting 
solutions  to  (19)  are  compared  to  the  exact  solution  for  a 
special  case. 


Table  1.  Height  adjustment  for  a four-hop  mode. 

Column  1 lists  the  pass,  or  break  point  index  (1  is 
Washington),  columns  2 and  3 list  north  latitude  and 
east  longitude  of  equally  spaced  break  points.  The 
fourth  column  shows  altitudes,  followed  by  break  point 
coordinates  adjusted  by  the  algorithm  to  account  for 
differing  hr.  The  seventh  column  lists  ground  angles 
which  correspond  to  the  adjusted  coordinates,  followed 
by  corresponding  pass  range  and  the  error  in  break  point 
position  due  to  our  approximation.  All  angles  are  in 
degrees  and  all  distances  in  km. 


1 Break 
Point 

Equal  Spacing 

r 

Adjusted  for 

hr 

Pass 

Range 

Position' 

Error 

Lat. 

Long. 

hr 

Lat. 

Long. 

Angle 

1 

33.90 

-77.03 

0 

38.90 

-77.03 

6.97 

1260.81 

0 

2 

44.84 

-65.19 

290.30 

45.20 

-64.28 

1260.81 

1 .88 

3 

49.30 

-51.01 

0 

49.77 

-48.86 

6.96 

1310.66 

3.76 

4 

51.74 

-34.72 

307.53 

51 .98 

-30.43 

1310.66 

6.43 

5 

51.75 

-17.54 

0 

51.15 

-11.44 

6.92 

1598.57 

9.10 

6 

49.35 

-1.23 

417.43 

j 46.34 

9.22 

1598.57 

6.59 

7 

44.91 

13.00 

0 

! 38.51 

25.65 

7.05 

548.21 

4.08 

8 

38.99 

24.86 

92.67 

j 35.37 

30.39 

548.21 

2.04 

9 

32.06 



34.76 
■ ■ ■ 

0 

; 32.06 

34.76 

0 
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The  mode  illustrated  has  4-hops,  three  F2-layer  reflections 
and  one  E-layer  reflection.  As  the  ray  traverses  the  great  circle  path 
successive  reflections  occur  at  increasing  effective  altitudes  until  the 
last  reflection,  which  is  at  sharply  decreased  altitude.  The  reasons  for 
this  behavior  illustrate  the  complicated  interplay  of  ionospheric  behavior 
and  radio  propagation.  The  second  reflection  point  is  above  the  first 
predominantly  because  the  F2  and  FI  layers  are  weaker  at  the  second  hop 

(break  point  4 in  Table  1)  where  the  local  time  is  about  4:15  a.m.  The 

third  hop  (break  point  6)  has  higher  effective  altitude  despite  the  dawn 
buildup  of  the  F2  layer,  because  the  FI  layer  is  nearly  strong  enough  to 
reflect  the  ray  and  thus  refracts  strongly.  Finally  at  the  fourth  hop 
(break  point  8)  not  only  is  the  FI  layer  capable  of  reflecting  the  ray 
but  also  the  E-layer  can  reflect  and  our  algorithm  chooses  the  lower  mode. 

An  exact  solution  would  have  yielded  equal  angles  approximately 
7.008°.  Thus  the  seventh  column  of  the  table  can  be  used  as  a measure  of 
inaccuracy  of  our  algorithm.  Another  measure,  one  which  is  more  directly 
pertinent  to  nuclear  effects,  is  the  error  in  placement  of  the  break  points. 
This  error  measures  the  distance  by  which  a disturbance,  beta  patch  or 

other,  might  be  missed  and  is  shown  in  the  final  column.  The  error  is 

seen  to  build  up  to  about  9 km  then  fall  back  toward  zero,  despite  the 
fact  that  break  points  were  shifted  as  much  as  1397  km  (No.  7).  Since 
any  natural  or  nuclear  disturbance  of  interest  will  be  large  compared  to 
9 km  one  concludes  that  the  approximation  is  adequate. 
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5. 


AMBIENT  IONOSPHERIC  TILT  ANGLE  CALCULATION 


The  main  motivation  for  inclusion  of  ambient  ionospheric  tilt 
effects  is  to  establish  a mechanism  by  which  tilt  effects  in  a highly  dis- 
turbed ionosphere  can  be  calculated.  For  this  reason  we  neglect  lower 
refracting  layer  effects  and  treat  only  reflecting  layer  tilt. 


In  the  sketch  the  solid 
lines  illustrate  real  ray  paths, 
one  of  which  connects  points  (1) 
and  (2) . This  ray  reaches  a max- 
imum height,  hm>  and  is  to  be 
represented  by  an  equivalent 
triangular  path  which  reflects 
at  height  hr>  illustrated  by 
dotted  lines.  To  find  the  effec- 
tive tilt  angle,  t,  of  the 
reflection  point,  we  calculate 


the  true  maximum  height,  h',  of  a similar  ray,  with  points  (1)  and  (2) 
displaced  a distance  6r. 


m 
Then 


T = 


(h  -m 

m m 
6r 


(30) 


Equation  30  is  good  to  first  order  which  is  adequate  for  the  first 
treatment  of  this  effect;  we  do  not  hope  for  accuracy  when  t is  large. 

For  the  moment,  and  perhaps  forever,  it  will  have  to  be  sufficient  to  know 
that  the  mode  is  highly  disturbed  when  T is  large,  and  the  approximate 
degree  of  disturbance. 


28 


5.1 


Tilt  Parallel  to  the  Great  Circle  Path 


After  calculation  of  the  first  approximation  to  reflection,  de- 
scribed in  Section  3,  the  true  maximum  ray  penetration  height,  h^,  is 
known  from  Equation  4 and  the  corresponding  value  of  x from  Equation  1. 

It  only  remains  to  choose  one  or  more  other  points  along  the  great  circle 
path  to  calculate  h^,  using  the  same  value  of  x.  Since  the  ionosphere 
is  described  by  parameters  specified  at  points  along  the  great  circle,  the 
technique  adopted  is  to  locate  the  two  great  circle  ionospheric  data  points 
surrounding  the  reflection  point,  evaluate  h^  at  these  two  points,  then 
calculate  "parallel  tilt",  T„  from  (30).  If  the  ray  is  capable  of  pene- 
trating the  ionosphere  at  one  of  the  great  circle  data  points,  the  reflec- 
tion point  is  used  in  its  place.  The  convention  is  adopted  that  tilt  is 
positive  if  height  increases  toward  the  receiver. 

5.2  Tilt  Perpendicular  to  the  Great  Circle  Path 

It  would  seem  that  our  method  of  ionospheric  description,  which 
restricts  knowledge  of  ionospheric  parameters  to  those  along  the  great 
circle  path,  might  prevent  finding  tilt  angles  perpendicular  to  the  great 
great  circle.  This  is  not  the  case  because  at  each  point  on  the  great  circle 
path  the  ionospheric  data  table  covers  all  24  hours  of  local  time.  This  is 
equivalent  to  data  at  all  values  of  longitude  provided  effects  of  geomag- 
netic coordinate  changes  on  the  longitude  dependence  in  the  neighborhood 
of  the  point  of  interest  are  neglected. 

The  situation  is  depicted  ^2 
in  the  sketch.  We  know  the  latitude, 
longitude  coordinate  of  the  reflec-  X 
tion  point  (X,Jl),  and  the  two  great 
circle  ionospheric  data  points 
surrounding  it  (A^,£^)  and  ^1 

^2 ' ^2 ) * t^lt  angle 

GC  Point  1 
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perpendicular  to  the  great  circle  path  can  be  obtained  by  differencing 

h at  the  points  (A,,£  ) and  (A„,5,  ). 
m r 1 e 2 p 

Since  arc  lengths  between  great  circle  data  points  are  small,  plane 
geometry  will  suffice.  In  this  approximation  the  interesting  longitudes  are 
given  by 


= l + (X-X1)CX2-A1)/(£2-i1) 


i = i * (X-X2)(A2-X1)/(£2-£1) 


(31) 

(32) 


If  the  local  time  at  (A,£)  is  t,  then  ionospheric  parameters 


for  (A^Jlg)  can  be  obtained  by  evaluating  those  at  (X ^ , £ ^ ) at 


t = t + a (£.-£) 


(33) 


where  a is  the  factor  appropriate  to  convert  longitude  to  time  (in  our 
case  ionospheric  time  is  measured  in  radians  from  local  midnight,  so  a = 1) 
Similarly 


t = t + a(Z  -Z) 
P P 


is  the  local  time  appropriate  for  point  (^2’^p-* 


(34) 


Then  the  perpendicular  tilt,  T L , is  obtained  from 


Ti  = thm(e)  ' hm(P)l/d 


(35) 


where  d is  the  distance  between  points  (X^,2.g)  and  (^2’^p^‘  Note  the 
convention  adopted  here  is  that  perpendicular  tilt  is  positive  if  height 
increases  to  the  right  when  looking  from  transmitter  toward  receiver. 
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There  are  two  special  cases.  If  the  path  segment  1 ■+  2 is  very 

nearly  North-South  < 10  rad)  we  evaluate  h^  at  (X,£  + 20 

minutes)  and  (A,£  - 20  minutes).  If  the  path  segment  1 -►  2 lies  too 

2 

nearly  East-West  (IX^-X^I  < 10  rad)  this  scheme  fails  and  zero  perpen- 
dicular tilt  is  assumed. 


6. 


EFFECT  OF  IONOSPHERIC  TILT  ALONG  THE  GREAT  CIRCLE  PATH 


Consider  the  case  of  a propagating  mode  with  arbitrary  but  fixed 
number  of  hops  where  the  ionosphere  varies  along  the  path  such  that  each 
hop  has  a different  virtual  reflection  height,  h^,  and  has  a small  tilt 
angle,  x^;  but  the  ray  remains  in  the  great  circle.  As  before,  the  path  is 
composed  of  straight  line  segments  which  are  called  "passes"  connected  at 
"break  points".  Thus  an  n-hop  mode  consists  of  2n  passes  and  2n  + 1 break 
points  as  depicted  in  the  sketch  below,  for  a 4-hop  mode. 


Break  Point  1 

© 

Transmitter 


Break  Point  9 

® 

Recei ver 


The  tilt  angles  are  given  by 

x - 

i dx  . 


(36) 


The  entire  mode  geometry  has  been  calculated  as  previously  de- 
scribed including  effects  of  varying  reflection  height,  but  under  the 
assumption  that  the  ionosphere  is  circularly  symmetric  about  the  center 
of  the  earth  at  each  reflection  point,  that  is,  effects  of  tilt  have  been 
neglected.  A formalism  is  needed  to  modify  break  point  coordinates  to  ac- 
count for  tilt,  one  which  is  accurate  for  reasonably  small  deflections  and 
"sensible"  for  large.  By  "sensible"  we  mean  the  break  point  range  should 
be  mom  tone  increasing  with  break  point,  reflection  altitude  shouldn't 
suffer  a major  fractional  change,  etc. 
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6.1 


Single  Hop  - Plane  Geometry 


Initially  consider  the  problem  of  a single  hop  in  plane  geometry, 
then  one  can  work  toward  a multihop  problem  in  circular  geometry.  In  the 
sketch  below,  dashed  lines  show  the  situation  which  would  obtain  if  the 
ionosphere  were  parallel  to  the  surface  of  the  earth  while  solid  lines 
show  the  situation  which  obtains  when  the  ionosphere  has  a (positive)  tilt 
angle  T. 


From  the  sketch,  the  launch  angle  without  tilt  is 

ctng  = tan  a = ^v-  (37) 

o 2h  ' 

and,  if  launched  at  angle  6 , the  ray  will  be  land  beyond  the  receiver  at 
angle  6-2x. 
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ctn  (3-2x) 


(38) 


As  a first  step  in  finding  the  adjustment  in  £,  h' , and  r 
required  to  force  the  ray  to  reach  the  receiver,  search  for  an  expression 
for  the  angle 


a = (Yj  + Y2)/2 


(39) 


as  a function  of  tilt  angle  x. 


Bisection  of  y by  the  dashed  line  shows 


y.  = a - x and  y~  = a + x 


(40) 


One  useful  relation  can  be  found  by  combining 


h'tan  (a-x)  = D/2  - 6 


(41) 


and 


h'tan  (a-x)  = D/2  + 6 

as 

h'[tan  (a+x)  + tan  (a-x)]  = D 
In  order  to  eliminate  h'  another  relation  is  required, 
tanx  = (h-h')/6 
with  (41)  to  eliminate  6, 


(42) 


(43) 

Use 

(44) 
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(h-h")  ctnx  = j - h'  tan  (a-x)  (45) 


Elimination  of  h between  (45)  and  (43)  followed  by  an  unusually  tedious 
amount  of  trigonometric  manipulation  leads  to  the  startling  result 

h tan  a = D/2  ! (46) 

In  short,  the  angle  by  which  the  successful  ray  is  reflected  by  the  tilted 
ionosphere  is  equal  to  the  untilted  value,  independent  of  X.  At  risk  of 
belaboring  the  obvious  we  point  out  that  this  result  is  not  a small  angle 
approximation  but  is  exact. 


Having  established  that  the  tilted  and  untilted  cases  have  a 
common  apex  angle  it  is  a short  exercise  in  plane  geometry  to  show  that 
the  launch  angle  has  been  raised 


SBt  = x 


(47) 


while  the  receiving  angle  has  been  decreased  by  x. 


From  (43),  now  that  a is  known  to  be  constant,  find  the  change 
in  reflection  altitude. 


2h  tana 


h'-  h 


tan (a  +x)  + tan (a  -x) 
o o 


(48) 


After  more  exercise  in  trigonometric  substitution,  Equation  48  simplifies  to 

. 2 


h - h = 


"m)‘ 


(49) 
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Using  (44)  one  finds  the  ground  range  of  the  reflection  point  to 
have  been  shifted  by  an  amount 


h-h' 

tanr 


sint  cost 
sin23 


(50) 


where  the  minus  sign  is  adopted  to  show  that  positive  X (reflection  alti- 
tude increasing  with  ground  range)  produces  a shift  toward  the  transmitter. 


This  single  hop  case  will  be  useful  primarily  as  a check  on  approx- 
imate methods  to  be  developed  for  more  complex  cases.  Attempts  at  preci- 
sion for  more  complex  cases  rapidly  become  bogged  down  in  more  trigonometry 
than  justified  on  the  basis  of  our  assumption  of  constant  x in  the 
vicinity  of  a given  break  point. 

6.2  Heights  Constant  - One  Place  Tilted  - Plane  Geometry 


Now  consider  the  cumulative  effects  of  multiple  hops  on  break 
point  range.  In  the  sketch  a ray  has  been  launched  at  angle  G appro- 
priate for  a multiple  hop  path  to  the  receiver  at  (R)  , the  nt*‘  break 
point . 
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The  tilt  T causes  reflection  angles  to  change  by  2t  after 
point  m.  Thus  the  distance,  6r,  by  which  the  receiver  is  missed  is 

6r  = (n-m)  h(ctn(B-2T)  - ctnB) 

**  (n-m)  h(l+ctn26)2T  . (51) 

The  total  range  from  (t)  to  (R)  is 

R = (n-1)  h ctnB  . (52) 


then 


HR  ? 

gg-  = - (n-l)h(l+ctn  B)  . (53) 

So  the  change  in  launch  angle  required  to  bring  the  ray  down  on  the  receiver 
is  found  from  requiring 

§ 66  - - 6r 


and  is 


dB 


2 (n-m) 
(n-1) 


(54) 


We  can  check  our  approximation  by  noting  that  for  a one-hop  case, 
n = 3,  m = 2,  (54)  gives 


dB  = t , the  exact  answer,  while  for  an  extreme  case,  m 

small , n 


dB  -*  2t, 
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Reference  to  the  sketch  shows  the  angle  2t  to  be  the  one  which  yields 
the  unperturbed  angle  3 thereafter.  For  this  extreme  case  the  approxi- 
mation has  introduced  an  error  equivalent  to  the  miss  distance  on  one  hop, 
an  error  commensurate  with  current  understanding  of  the  ionosphere,  especi- 
ally so  because  some  of  the  error  can  be  eliminated  by  retroactive  adjust- 
ment of  break  points  to  force  the  ray  to  land  on  the  receiver. 


The  increment  in  range  at  intermediate  break  points  is  given  by 


h ctn  [3  + 


2 Tl 
(n-1)  TJ 


i ^ m 


fir . = 
1 


h ctn  [3 


2 0l-.1I  t1 

(n-1)  TJ 


i > m 


(55) 


6.3 


Heights  Vary  - One  Place  Tilted  - Plane  Geometry 


To  find  the  miss  distance  we  must  sum  heights  in  Equation  51 


fir  = T.  h.  (l+ctn^B)  2t 


2 (1+ctn  B)t  £ h -h 

m+1  1 1'1 


(56) 


Following  an  argument  similar  to  the  previous  case  we  find 
n 


63  = t 


2 V h.-h. 
m+1  1 1'1 


y Ih.-h.  . 
2 I 1 1-1 


(57) 
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6.4 


Heights  Vary  - All  Places  Tilted  - Plane  Geometry 


From  (57)  the  change  in  launch  angle  due  to  tilt  at  location 


m is 


60  = T 

m m 


n 

2 £ hi 
i=m+l 
n 

£ \ 

i=2 


(58) 


The  total  change  in  launch  angle  is  given  by  summing  (58) ; 

n 


60  = 


n-1 

2 £ 2>i 

m=2  i=m+l 


n 

m=  2 1 


(59) 


The  simplest  method  of  finding  the  new  ranges  and  heights  of 
break  points  is  to  accumulate  them,  viz. 
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(6 


Sh  = (r  -r  )tan  x 
m mm  m 


where  a second  order  effect  in  the  equation  for  r has  been  neglected. 

m 


When  one  reaches  the  last  break  point  the  accumulated  error  is 
measured  by  the  miss  distance.  All  ranges  are  then  scaled  to  force  agree 
ment  in  total  range. 

6.5  Extension  to  Circular  Geometry 

If  3 is  interpreted  as  the  launch  angle  relative  to  the  chord 
of  the  arc  as  depicted  below  nothing  changes  to  first  order  in  the  pass 
range  over  the  radius  of  the  earth. 


It  is  necessary,  however,  to  rederive  (59)  accounting  for  the  fact 

that,  while  incident  and  reflection  angles  are  equal  relative  to  the  earth's 

surface,  they  are  not  equal  relative  to  the  chord  of  the  arc.  The  effect  of 

2 

this  extension  is  that  the  factor  (1+ctn  3)  cannot  be  regarded  as  constant, 
but  must  be  retained.  Equation  59  becomes 


63 


n-1 

L 

m=2 


m 


i=m+l 


Ik  (1+ctn  3^) 


L h.(l+ctn23i) 
i=2 


(62) 


In  many  cases  a second  order  term  due  to  the  change  in  incident 
angle  which  results  from  change  in  pass  range  is  significant.  Proper  treat- 
ment of  this  term  is  a bit  cumbersome  and  more  importantly,  tends  to  de- 
stabilize the  adopted  mathematical  treatment.  A satisfactory  expedient  is 
to  retroactively  adjust  all  pass  ranges,  to  force  agreement  with  total  range, 
by  an  amount  proportional  to  the  square  of  the  calculated  pass  range. 


6.6  Example 


Table  2 summarizes  results  of  the  approximate  treatment  developed 
above  when  applied  to  the  same  four-hop  Washington  to  Tel  Aviv  mode  illustrated 
previously.  This  particular  mode  was  chosen  primarily  because  it  represents 
as  difficult  a case  for  the  parallel  tilt  approximation  as  has  been  encountered 
to  date.  The  table  shows  errors  as  large  as  about  20  km  in  placement  of 
break  points.  Since  this  error  is  small  compared  to  the  size  of  any  possible 
D-layer  absorbing  region  the  method  is  considered  satisfatory.  The  calcu- 
lated launch  angle  at  Washington  is  5.11°  compared  to  the  exact  solution  of 
5.01°  and  the  calculated  received  angle  at  Tel  Aviv  is  12.15°  compared  to 
the  exact  solution  of  11.89°. 
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Table  2.  Parallel  tilt  adjustment  for  a 4-hop  mode. 


Column  1 lists  the  break  point  index  (1  is 
Washington,  D.C.,  9 is  Tel  Aviv).  The  next 
three  columns  list  latitude,  longitude,  and 
great  circle  range  measured  from  Washington 
after  adjustment  for  different  reflection 
height.  The  following  four  columns  show  the 
tilt  angle  parallel  to  the  great  circle  path, 
reflection  height  and  position  adjustments 
for  tilt  and  the  increase  in  range  due  to  tilt. 
The  next  column  shows  the  exact  trigometric 
position  (accepting  h~)  and  the  last  column 
the  amount  by  which  tne  approximate  method 
overestimates  the  position.  All  distances 
are  in  km  and  all  angles  in  degrees. 


7.  TILT  ACROSS  GREAT  CIRCLE 

7.1  Plane  Case 

Looking  from  transmitter  toward  receiver,  one  sees  a tilted  ionos- 
phere and  the  ray  path.  Positive  tilt  means  h increases  to  the  right  as 
depicted  below.  Launch  a ray  along  the  great  circle  path.  Without  tilt 

one  would  see  only  one  vertical  line.  With  tilt  n at  break  point  m 

m r 

the  ray  is  deflected  to  the  right. 


The  ray  comes  off  the  ionosphere  at  angle  2n  , causing  it  to  land  a distance 

m 

21^  to  the  right  of  the  great  circle  path,  then  if  there  are  no  more 

tilted  break  points,  it  moves  additional  distance  2p  h.  for  each  pass. 

m j r 


The  total  displacement  at  the  end  of  the  path,  break  point  n, 
due  to  the  tilt  at  break  point  m is 


d 

m 


2n 


m 


n 


£ |h.-h 

i=m+l ' 


(63) 


If  one  assumes  a tilt  at  each  break  point,  one  has 
n-1  n-1  n 

d = Z d = 2 £ 1 V lli.-h.  ,1 

t m m . *-•  . l l-l 

m=2  m=2  i=m+l‘ 
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The  launch  azimuth  (positive  to  right)  of  the  ray  which  reaches 
the  receiver  must  counteract  the  deflection; 


aL  - - d/R 


where  R is  the  great  circle  range  of  the  total  path. 

The  azimuth  of  the  ray  at  the  receiver  is 


a = - d . / r 
R n-1  n 


where  r^  is  the  ground  range  of  the  final  pass. 


Since  the  launch  azimuth  is  a. , the  resulting  displacement 
perpendicular  to  the  great  circle  of  the  J break  point  is 


J-l 

J 

sT  = 2 Y n 

Y h.-h.  + ad 

.*->  l l-l  L J 

(67) 

m=2 

i=m+l  1 1 

where  d.  is  ground  range  to 

point  J. 

The  new  reflection  altitude  is 


hJ  = hJ  + SJ  tan  nj 


n J-l 

Numerically  it  is  convenient  to  collect  Y |h.-h.  ,1  and  Y P , then 

i=m+l  1 1 m=2 

s,  * sj.i*  ">j-i  ?,  ihi-hi-ii  * 2ihJ-hJ-iiMn»  * “l  rj 


i -h  T , V n + a.  r 
J J-l  m L J 


where  r^  is  the  incremental  range  of  the  J pass. 


7.2  Extension  to  Spherical  Geometry 

If  Ik  is  interpreted  as  the  height  of  reflection  above  the  chord 
of  the  hop  rather  than  above  the  earth's  surface.  Equation  69  is  still  valid. 

7.3  Determination  of  Break  Point  Horizontal  Coordinates 


If  A is  colatitude  and  £ is  longitude,  to  a good  approxima- 
tion one  can  write 


Sj  = Re  ^(<5A)2  + sin2  Aj  (6£)‘ 


(70) 


Since  the  slope  of  the  line  segment  Sj  is  the  negative  gradient 


of  the  line  segment  rT, 

J 


dA 

d£ 


■ - (§)  - - 
'T 


-1  'V'j-i) 


(71) 


combining  (70)  and  (71), 

6A  = 


± S, 


(72) 
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and. 


where  the  positive  sign  is  chosen  if  £ > £ , 

J “ I 


J * sin2xj  t73) 

where  the  positive  sign  is  chosen  if  Aj  < Aj  . 


7.4  Example 

Table  3 illustrates  results  obtained  for  the  Washington-Tel  Aviv 
case.  The  distance  by  which  the  ray  missed  the  receiver  (1.5  parts  in  10  ) 
measures  the  combination  of  second  order  error  in  the  technique  and  round 
off  error  in  a 32  bit  computer  word  length.  Of  more  interest  is  the  much 
lower  magnitude  of  perpendicular  displacements  relative  to  parallel  dis- 
placements. Most  of  this  reduction  in  effect  is  due  to  the  circumstance 
that  the  path  lies  mostly  East-West  so  that  perpendicular  gradients  are 
smaller  than  parallel  gradients  and  the  further  accident  that  perpendicular 
tilts  tend  to  cancel  (the  first  two  angles  sum  to  +0.26  and  the  next  is 
-0.26).  However,  parallel  gradients  systematically  tend  to  be  about  an 
order  of  magnitude  more  effective  in  displacement  than  perpendicular  gra- 
dients. The  parallel  effect  tends  to  be  enhanced  by  small  angle  of  incidence 
at  the  earth,  while  the  perpendicular  effect  always  acts  as  though  the  ray 
were  near  a right  incident  angle. 
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Table  3.  Perpendicular  tilt  adjustment  for  a 4-hop  Mode. 

Column  1 lists  the  mode  break  point  index,  column  2 
reflection  altitude,  and  column  3 the  tilt  angle 
perpendicular  to  the  great  circle  path.  Column  4 
lists  displacement  of  break  point,  positive  to  right 
of  great  circle  path  as  viewed  from  the  transmitter, 
and  column  5 the  adjusted  reflection  height  due  to 
perpendicular  tilt.  All  distances  are  in  kilometers 
and  all  angles  in  degrees. 


Break 

Point 

1 


hr 

n 

After  Adjustment 

s 

h 

0 

0 

0 

0 

290.03 

0.11 

-1 .29 

290.03 

0 

0 

-1 .40 

0 

308.50 

0.15 

-1 .49 

308.49 

0 

0 

0.02 

0 

408.41 

-0.26 

2.34 

408.40 

0 

0 

0.91 

0 

92.54 

-0.03 

0.43 

92.54 

0 

0 

-0.14 

0 
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Code  7700,  Timothy  P.  Coffey 
Code  5410,  John  Davis 


Commander 

Naval  Space  Surveillance  System 
ATTN:  CAPT  J.  H.  Burton 

Naval  Space  System  Activity 
ATTN:  A.  B.  Hazzard 


Air  Force  Avionics  Laboratory,  AFSC 
ATTN:  AAB , H.  M.  Hartman 
ATTN:  AAD,  Wade  Hunt 
ATTN:  AAD,  Allen  Johnson 

Headquarters 

Electronics  Systems  Division/XR 
ATTN:  XRC,  Lt  Col  J.  Morin 
ATTN:  XRE,  Lt  Michaels 

Headquarters 

Electronic  Systems  Division/YS 
ATTN:  YSEV 

Headquarters 

Electronic  Systems  Division  (AFSC) 
ATTN:  Jim  Deas 
ATTN:  James  Whelan 
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DEPARTMENT  OF  ENERGY  (Continued) 


Commander 

Foreign  Technology  Division,  AFSC 
ATTN:  ETD,  B.  L.  Ballard 
ATTN:  NICD,  Library 


Sandia  laboratories 
Livermore  Laboratory 

ATTN:  Doc.  Con.  for  Thomas  B.  Cook 
ATTN:  Doc.  Con.  for  Byron  E.  Murphey 


Hq.  USAF/RD 

ATTN:  RDQ 


OTHER  GOVERNMENT  AGENCIES 


Commander 

Rome  Air  Development  Center,  AFSC 
ATTN:  EMTLD,  Doc.  Library 
ATTN:  V.  Coyne,  OCSE 

Commander 

Rome  Air  Development  Center,  AFSC 
ATTN:  ETEI,  A.  Lorentzen 

SAMSO/MN 

ATTN:  MNNL,  Lt  Col  Kennedy 


Department  of  Commerce 
Office  of  Telecommunications 
Institute  for  Telecom  Science 
ATTN:  G.  Reed 
ATTN:  L.  A.  Berry 
ATTN:  William  F.  Utlaut 
ATTN:  A.  Glenn  Jean 
ATTN:  D.  D.  Crombie 

NASA 

Goddard  Space  Flight  Center 

ATTN:  ATS-6  Ofc.,  P.  Corrigan 


SAMSO/SK 

ATTN:  SKA,  Lt  Maria  A.  Clavin 
SAMSO/SZ 

ATTN:  SZJ,  Maj  Lawrence  Doan 
ATTN:  SZ 


SAMSO/YA 

ATTN: 


Commander  in 


YAT,  Capt  L.  Blackwelder 
Chief 


Strategic  Air  Command 

ATTN:  XPFS,  Maj  Brian  G.  Stephan 
ATTN:  NRT 

ATTN:  ADWATE,  Capt  Bruce  Bauer 
ATTN:  Chief  Scientist 


National  Oceanic  & Atmospheric  Admin. 

Environmental  Research  Laboratories 
ATTN:  C.  L.  Rufenach 
ATTN:  Joseph  H.  Pope 
ATTN:  Richard  Grubb 

Department  of  Commerce 

National  Bureau  of  Standards 
ATTN:  Raymond  T.  Moore 

Department  of  Transportation 

Office  of  the  Secretary 
ATTN:  R.  L.  Lewis 
ATTN:  R.  H.  Doherty 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 


DEPARTMENT  OF  ENERGY 

University  of  California 

Lawrence  Livermore  Laboratory 

ATTN:  Tech.  Info.  Dept.  L-3 
ATTN:  Frederick  D.  Seward,  L-46 
ATTN:  Glenn  C.  Werth,  L-216 

Los  Alamos  Scientific  Laboratory 

ATTN:  Doc.  Con.  for  R.  F.  Taschek 
ATTN:  Doc.  Con.  for  Donald  R.  Westervelt 
ATTN:  Doc.  Con.  for  P.  W.  Keaton 

Sandia  Laboratories 


Aeronomy  Corporation 

ATTN:  S.  A.  Bowhill 

Aerospace  Corporation 

ATTN:  S.  P.  Bower 
ATTN:  V.  Josephson 
ATTN:  D.  P.  Olsen,  120  Rm.  2224E 
ATTN:  Irving  M.  Garfunkel 
ATTN:  Norman  D.  Stockwell 
ATTN:  F.  E.  Bond,  A1  Rm.  5003 
ATTN:  F.  A.  Morse,  A6  Rm.  2407 
ATTN:  T.  M.  Salmi 
ATTN:  J.  E.  Carter,  120  Rm.  2209 
ATTN:  SMFA  for  PWW 


ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


Doc.  Con.  for  W.  D.  Brown,  Org.  1353 

Doc.  Con.  for  A.  Dean  Thornbrough,  Org.  1245  Analytical  Systems  Engineering  Corp. 


Doc.  Con.  for  D.  A.  Dahlgren,  org.  1722 
Doc.  Con.  for  J.  P.  Martin,  Org.  1732 
Doc.  Con.  for  Space  Proj.  Div. 


ATTN:  Radio  Sciences 


ATTN:  Doc.  Con. 

for  3141,  Sandia  Rpt.  Coll. 

ATTN: 

Glen  Keister 

ATTN: 

D.  Murray 

Department  of  Energy 

ATTN: 

J.  F.  Kenny 

Albuquerque  Operations 

Office 

ATTN: 

Glenn  A.  Hall 

ATTN:  Doc.  Con.  for  D.  W.  Sherwood 


Department  of  Energy 

Division  of  Headquarters  Services 

ATTN:  Doc.  Con.  for  Allen  Labowitz 


Brown  Engineering  Company,  Inc. 

ATTN:  Romeo  A.  Deliberis 

University  of  California  at  San  Diego 
Marine  Physical  Lab.  of  the  Scripps 
Institute  of  Oceanography 

ATTN:  Henry  G.  Booker 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Contir 

Charles  Stark  Draper  Laboratory,  Inc. 

JAYCOR 

ATTN 

J.  P.  Gilmore,  MS  63 

ATTN 

S.  R.  Goldman 

ATTN 

D.  B.  Cox 

Johns  Hopkins  University 

Computer  Sciences  Corporation 

Applied  Physics  Laboratory 

ATTN 

H.  Blank 

ATTN 

Thomas  Potemra 

ATTN 

John  Dassoulas 

COMSAT  Laboratories 

ATTN 

Document  Librarian 

ATTN 

R.  R.  Taur 

Linkabit  Corporation 

Cornell  University 

ATTN 

Irwin  Jacobs 

Department  of  Electrical  Engineering 

ATTN 

D.  T.  Farley,  Jr. 

Lockheed  Missiles  & Space  Co.,  Inc. 

ATTN 

Dept.  60-12 

ESL,  Inc. 

ATTN 

D.  R.  Churchill 

ATTN 

James  Marshall 

ATTN 

J.  Roberts 

Lockheed  Missiles  & Space  Co.,  Inc. 

ATTN 

C.  W.  Prettie 

ATTN 

Martin  Walt,  Dept.  52-10 

ATTN 

V.  L.  Mower 

ATTN 

Richard  G.  Johnson,  Dept.  52 

ATTN 

W.  L.  Imhof,  Dept.  52-12 

Ford  Aerospace  & Communications  Corp. 

ATTN 

J.  T.  Mattingley,  MS  X22 

MIT  Lincoln  Laboratory 

ATTN 

James  H.  Pannell,  L-246 

General  Electric  Company 

ATTN 

D.  Clark 

Space  Division 

ATTN 

Lib.  A-082  for  David  M.  Towl 

Valley  Forge  Space  Center 

ATTN 

Mr.  Walden,  X113 

ATTN:  M.  H.  Bortner,  Space  Sci.  Lab. 

McDonnell  Douglas  Corporation 

General  Electric  Company 

ATTN 

Wi 1 1 iam  01  son 

TEMPO-Center  for  Advanced  Studies 

ATTN 

N.  Harris 

ATTN 

Tim  Stephens 

ATTN 

J.  Moule 

ATTN 

Mack  Stanton 

ATTN 

George  Mroz 

ATTN 

Warren  S.  Knapp 

ATTN 

B.  Gambill 

University  of  Lowell  Rsch.  Foundation 

ATTN 

DAS I AC 

ATTN 

Dr.  Bib! 

ATTN 

Don  Chandler 

Mission  Research  Corporation 

General  Electric  Company 

ATTN 

Dave  Sowle 

ATTN 

F.  A.  Reibert 

ATTN 

F.  Fajen 

ATTN 

R.  Bogusch 

General  Research  Corporation 

ATTN 

R.  Hendrick 

ATTN 

John  Ise,  Jr. 

ATTN 

P.  Fischer 

ATTN 

Joel  Garbarino 

ATTN 

M.  Scheibe 

ATTN 

Steven  L.  Gutsche 

Geophysical  Institute 

5 cy  ATTN 

Tech.  Lib. 

University  of  Alaska 

ATTN 

Neal  Brown 

The  Mitre  Corporation 

ATTN 

Technical  Library 

ATTN 

Chief  Scientist,  W.  Sen 

ATTN 

T.  N.  Davis 

ATTN 

C.  E.  Callahan 

ATTN 

G.  Harding 

GTE  Syl vania,  Inc. 

ATTN 

John  Moganstern 

ATTN 

Marshal  Cross 

Pacific-Sierra  Research  Corp. 

HSS,  Inc. 

ATTN 

E.  C.  Field,  Jr. 

ATTN 

Donald  Hansen 

Pennsylvania  State  University 

University  of  Illinois 

Ionosphere  Research  Laboratory 

Department  of  Electrical  Engineering 

ATTN 

Ionospheric  Research  Labs. 

ATTN 

K.  C.  Yeh 

Photometries,  Inc. 

Institute  for  Defense  Analyses 

ATTN 

Irving  L.  Kofsky 

ATTN 

Joel  Bengston 

ATTN 

J.  M.  Aein 

Physical  Dynamics,  Inc. 

ATTN 

Ernest  Bauer 

ATTN 

E.  J.  Fremouw 

ATTN 

Hans  Wolfhard 

Physical  Dynamics,  Inc. 

Inti.  Tel. 

& Telegraph  Corporation 

ATTN 

A.  Thompson 

ATTN:  Tech.  Library 

ATTN 

Joseph  B.  Workman 
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R&D  Associates 

ATTN:  Bryan  Gabbard 
ATTN:  Robert  E.  Lelevier 
ATTN:  William  J.  Karzas 
ATTN:  William  B.  Wright,  Jr. 

ATTN:  Forrest  Gilmore 
ATTN:  H.  A.  Ory 
ATTN:  R.  P.  Turco 

The  Rand  Corporation 

ATTN:  Cullen  Crain 
ATTN:  Ed  Bedrozian 

Raytheon  Company 

ATTN:  Barbara  Adams 

Science  Applications,  Inc. 

ATTN:  Jack  McDougall 
ATTN:  Lewis  M.  Linson 
ATTN:  D.  Sachs 
ATTN:  Curtis  A.  Smith 
ATTN:  E.  A.  Straker 
ATTN:  Daniel  A.  Hamlin 

Science  Appl ications,  Inc. 

Huntsville  Division 

ATTN:  Dale  H.  Divis 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

SRI  International 

ATTN:  M.  Baron 
ATTN:  L.  L.  Cobb 
ATTN:  Ray  L.  Leadabrand 
ATTN:  Walter  G.  Chesnut 
ATTN:  G.  Smith 
ATTN:  Walter  Jaye 
ATTN:  Alan  Burns 
ATTN:  Donald  Neilson 
ATTN:  David  A.  Johnson 
ATTN:  Charles  L.  Rino 
ATTN:  Gary  Price 
ATTN:  George  Carpenter 
ATTN:  James  R.  Peterson 

System  Development  Corporation 
ATTN:  E.  G.  Meyer 

TRI-COM,  Inc. 

ATTN:  Darrel  Murray 

TRW  Defense  & Space  Sys.  Group 

ATTN:  R.  K.  Plebuch,  Rl-2078 
ATTN:  Saul  Altschuler 
ATTN:  Diana  Dee 

Vi  si  Dyne,  Inc. 

ATTN:  J.  W.  Carpenter 
ATTN:  Charles  Humphrey 
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